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ABSTRACT 


The article presents the theoretical studies results of 
hydraulic shock with flow continuity rupture in the 
pressure the pumping station pipeline. Calculated 
dependencies are obtained for calculating the shock 
pressure magnitude and the first period duration of 
reduced pressure, taking into account the hydraulic 
resistance along the pipeline length in the discontinuity 
event in the flow. The calculation results are compared with 
the Professor D.N. Smirnov’s experiments results. 


Comparative calculations using the proposed method with 
D.N. Smirnov experimental data gives a satisfactory 
agreement and confirms the proposed method reliability. 
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1. INTRODUCTION 

Currently, in the Republic of Uzbekistan, 1,693 pieces. of 
irrigation pumping stations are used to irrigate more than 
2.3 million land hectares. Ensuring the normal operation 
safety and reliability of the pumping stations used is in 
current importance. A dangerous process during normal 
operation of pumping stations pressure pipelines, as you 
know, are hydraulic shocks, accompanied by the breaks 
formation in the flow continuity in the injection water 
supply. The water hammer study, taking into account the 
continuous flow rupture in pressure systems, has not yet 
been fully studied [1,2,3,4]. Any scientific information on the 
water hammer study is relevant. This work is devoted to this 
issue. 


An experimental water hammer study taking into account 
the flow discontinuity is given in [3,4,5,6,7]. These works 
[3,4,5,6,7] did not take into account the reduced pressure 
duration in the first impact phase, no dependencesare given 
to determine the phases number of the water hammer. This 
data is very important for the impact parameters analysis. 


In [8], theoretical studies have been carried out, which are 
related to the discontinuity formation in the water flow ina 
pipeline and by determining the maximum pressure 
magnitude during hydraulic shock in different sections along 
its length. However, in [8], short pipelines are provided, in 
which the pressure losses due to friction along the length can 
be ignored. 


The work purpose is to determine the maximum pressure 
during a hydraulic shock in the first period of reduced 


pressure in the flow rupture event continuity, taking into 
account the pressure loss due to friction in a long pressure 
pipeline. 


2. CALCULATION METHOD 

Water is supplied by a pump at a speed Vo to the pressure 
basin. A non-return valve is installed at the beginning of the 
discharge pipeline. When the check valve is instantly closed, 
the water stops and the water pressure instantly decreases 
by the Vo value, where g is the gravity acceleration, ais the 

& 
propagation velocity of the shock wave. 


With the above hydraulic process, the following condition is 
met: 


OE oe +H (1) 


6aK.MQKC. mp 


where Heax. max. IS the value of the vacuum in the source of the 
hydraulic shock; 


Hrp. is friction head loss up to the considered point of the 
pressure pipeline cross-sections; 


Ho is geodetic head. 


When condition (1) is met, the flow continuity rupture 
occurs at the check valve, and water continues to move ata 
certain residual velocity through the pressure pipeline. 


The Allievi chain equations connecting the heads Hj and H;-1 
- and the velocities V; and V;-1 - in the section of the pipeline 
at the gate in the i-th and (i-1) - th phases have the form 
[9,10,11]: 


a 
H,_,+H,-2H,=—VV,-V,,), (2) 
§ 
for the first phase Hi = - Hsax. max. 


Then from (2) we obtain: 
A,-HA -2H,=—(V,-\), 
& 


6BaK.MQKC. 


from where 


ence ), Vi -represents 


V1 =Vo-V+where V, = S(H, +H 
a 
the residual speed with which the water mass torn off near 
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the check valve moves in the pipeline towards the pool. This 
speed decreases from phase to phase (friction forces are not 
taken into account here) until it reaches zero. 


The check valve has a vacuum space filled with air bubbles 
and steam, the pressure in it will be equal to, ie. H = - Hpax. 
max, the velocity of the separated water mass of the flow in 
the subsequent phases according to (2) will accordingly be 
equal to: 


V,=V,-V. 
V, =V, -3V. 
V, =V) —5V, (4) 


V, =V, —(2i-1V. 


We determine the total values of the traversed paths 
of the water mass for each phase, taking into account that 
this path first increases to a certain value S, and then 
decreases to zero. Then, equating to zero the indicated sum, 
we find the number of vacuum spaces 0, i.e. 


Wy, =0 
i=1 


or 
DV, =n -nV.)=0, (5) 
i=l 
From it 
Vo 
pe. 6 
V. (6) 


Substituting the value of n in (4), we obtain the velocity of 
the water mass at the moment of its impact on the check 
valve, equal to 

Vn = -(Vo - V4) = -V1. (7) 


In this case, the additional pressure of the water hammer in 
excess of the geodesic pressure will be equal to 


i ae ok SS, (8) 
& 


where AH; = Ho + Heax.max. the largest drop in shock pressure 
compared to geodesic; 


AF 406. =“y, =“(v,-¥.). (9) 
&§ & 


The largest value of the shock pressure will be equal to: 


H an tam, ao. (10) 
§ 


MaKc. 


When the water supply through the pump is instantaneously 
interrupted with the formation ofa discontinuity in the flow 
at the check valve in the pipeline, an increased pressure 
above the geodetic pressure arises, equal to [8]: 


AH, = 20, 


& 


The duration of the first period of reduced pressure without 
regard to friction forces is 
tv=n tg. (11) 


It is possible to determine over the geodetic pressure and 
the duration of the first period of reduced pressure, taking 
into account the head loss due to hydraulic resistance along 
the length of the discharge pipeline. We use the Darcy- 
Weisbach formula to determine the head loss 


pV re 
mp d 2¢ » 
where A - is the coefficient of hydraulic friction resistance, 
d - is the discharge line diameter. 


Then the loss of kinetic energy will be equal to: 
Lbs 
V,, = =H ay, 


mp mp 
a 


(13) 


At the moment of rupture of the continuity of the flow 
through the pipeline from the place of the rupture of 
propagation, the wave of reduced pressure (Ho + Haaxmax), 
which corresponds to the instantaneous velocity V:. 


This speed, formed due to the deformation of the walls of the 
pipeline and changes in the density of water in it, propagates 
through the pipeline. 


Calculating the loss in speed by formula (13), through each 
half-phase we will have 
V* = Vs - Vip. (14) 


where 


Vo =~ Ve. 
2a 


Obviously, for the entire time that there is a vacuum space 
(in the case under consideration, at the beginning of the 
pressure pipeline, at the check valve), there will be H= - 
Hysaxmax, and in the last section of the pressure pipeline 
adjacent to the pressure basin of sufficiently large 
dimensions positive pressure equal to atmospheric pressure 
plus the immersion of the last section under the pressure 
basin horizon. Based on this, in the initial section of the 
pressure pipeline (at the check valve), the instantaneous 
velocity during the rupture of the pressure flow changes by 
the value V:, and in the last section by the value V" = V«- V"rp 
(taking into account the speed loss for half phase of water 
hammer). 


Therefore, at the moment of restoring the continuity of the 
flow, the super geodetic pressure of the hydraulic shock will 
be equal to 


AH, =—V°+4Vv, =“(v'+v,), (15) 
g 8" Bg 


a ‘ 
where — V. the value of the pressure increase, equal to the 
& 


previous pressure decrease relative to the geodetic, taking 
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* a 
into account the loss inspeed V for half the phase; — V, is 


the value of the additional shock pressure obtained as a 
result of a water hammer a mass of water at a speed of 


V_ about a check valve. 


In this case, the maximum pressure of the water hammer 
will be equal 
a * 
H ge = Hy + AH, =H, +<V’ +V,). (16) 
& 


MaKc 


To determine the duration of the first period of reduced 
pressure and the velocity Vn of the water mass at the 
moment it hits the check valve, we will write expressions for 
the column velocity for each subsequent half of the water 
hammer phase. 


Residual speed value, i.e. the velocity of the mass of water at 
the beginning of the first phase is equal to: 
Vi = Vo = Vs, 


For the next half phase, this speed will be: 

Vi3 = Vi —- V+ - Vizp = Vo - 2V*- Virp, 

and at the beginning of the second phase, the speed value 
will be: 

V2 Vis - Ve - Vis Tp Vo - 3V«- (Vi tp + Vis tp) 

and so on. 


Thus, the velocities of the mass in the next half of the phase 
will respectively be equal: 


V,=V)-V., 
Vis = Vo ~~ 2V. = Vis > 
Vy = Vo =3V5 ~ (V,.,, +V, snp) 


, (17) 


V, =Vo — (20-1 We — (Vii + Vj 5p Tevet 


n lmp 1,5mp 


+Vin_o.s) Vamp = seals V(n-0,5)mp 


mp i 





here k - the number of the phase from which the mass of 
water begins its reverse movement, and the speed calculated 
by formula Virp (i = 1; 1,5; ...; n - 0,5) (13). 


Similarly to formula (5), the algebraic sum of all velocities 
(17) for the first period ofreduced pressure is equal to zero: 
(2n - 1)Vo - n(2n - 1)V* - (2n - 2)Virp - (2n - 3)Vi sep - (2n - 
A)Vorp - 

a [2n = 2(k = 0,5)]V¢x-0,5)rp + (2n = 2k) Vx tp + [2n = 2(k + 
0,5)]Vx+0,5)rp +... = 0 


or 
2V= n2 - (2Vo + V+- Q)n + (Vo - R) = 0, (18) 


where 

Q = 2(Virp + Visep oF Vik-0,5)7p = Virp i Vic+0,5)rp = i) 

R= 2Virp + 3V15rp +4 Varp Fait 2(k - 0,5) Vex-0,5)rp -2k Vierp 7 
- 2(k + 0,5)Veceo,s)rp - 


Solving the quadratic equation (18) with respect to n, we 
find the number of phases of a water hammer (during which 
there was a cavitation-vacuum space) equal to 


| 2 
AT ae +V, —Q)+¥(2V, +V.-O) -8V.(V -R)| 


(19) 

here only a positive value of the root is taken by equation 
(18), since the plus sign satisfies the condition of the 
problem being solved. 


Having found n, from (17) we can determine the water mass 
velocity at the moment of its impact on the check valve. 


3. STUDIES RESULTS 

Below are calculations to determine the magnitude of the 
shock pressure and the duration of the first period of 
reduced pressure, taking into account the hydraulic 
resistance along the length of the pipeline in the event of a 
rupture of the flow continuity. The calculation results are 
compared with the results of the experiments of D.N. 
Smirnov [6]. 


As can be seen from Table 1, comparative calculations by the 
proposed method with the experimental data of D.N. 
Smirnov give satisfactory results. 



















































































Table 1 
Ys 7 
110 : F 
Ne of according to according to 
F AH; No onbiTamM pacue- ; A 

SSH w TONG) iv NHN bw Cn aks Litre 

: CmMupHoro = TaM : ? 
experiments experiments calculations 

aBTopa 
M Nu 

1 0,42|36,5| 7,0 | - |1148| 39,4 | 922 |35,0 33,0 33,8 |2,49 2,85 2,69 

2 1,16 |40,0|55,0| 7 |1148/109,0; 922 |47,0 86,0 89,3 |2,50 5,10 5,10 

3 1,32 |21,0|70,5| 7 |1180/105,6| 787 |28,0 82,0 74,6 |3,00 8,10 8,22 

4 1,44 |17,0|75,0| 7 {1176} 91,2 | 622 |24,0 56,0 61,5 {3,80 9,00 9,72 

5 1,50 |10,0/82,5| 8 {1170} 83,5 | 545 |18,0 49,0 48,9  |4,80 11,10 13,92 

6 1,69| 5,0 |87,0| 8 {1165} 71,7 | 416 |13,0 32,0 36,0 [5,60 14,70 13,10 

7 1,22 |36,5|56,0| 7 |1148/114,7) 922 |43,5 88,5 89,3 |2,49 5,40 5,09 

8 0,93 |36,5/36,5| 6,5 |1148| 87,8 | 922 |43,0 64,5 71,4 {2,49 4,50 4,68 

9 0,76|37,0|22,0| 6 |1148| 71,3 | 922 |43,0 61,0 62,0  |2,49 3,90 3,96 

10 0,58 |37,0/14,0| 5 |1148| 54,3 | 922 |42,0 46,0 46,0  |2,49 3,45 3,19 

11 0,42 |37,0] 6,5 1148] 39,4 | 922 |35,0 33,8 2,49  |2,49 2,85 2,69 
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1. For experience II we have data: 
Vo = 1,16 m/sec; Ho = 40,0 m; Hrp = 55,0 M; Heacmax. = 7,0 m; a 
= 922 m/sec; ty = 2,50 sec; d= 50 mm. 


According to the initial data of the second experiment, we 
calculate: 


1 L 
aVo =109,0; =“ -—=0,44 sec/m; V+ = 0,50 m/sec; Vrp 


& a 
= 0,44-0,502 = 0,11 m/sec; 


V* = V«- V'trp = 0,50 - 0,11 = 39 m/sec; Vi = Vo - V+ = 1,16 - 
0,50 = 0,66 m/sec; 


Virp = 0,44-0,662 = 0,19 m/sec; Vi5 = Vi - V« - Virp = 0,66 - 
0,50 - 0,19 = 0,03 m/sec; 


Vi5tp = 0,44-0,032 = 0,00; V2 = Vi5 - V+ + Visrp = 0,03 - 0,50 + 
0,00 = 0,53 m/sec; 


Vorp = 0,44-0,532 = 0,12 m/sec; V25 = V2 - V+ + Varp = - 0,53 - 
0,50 + 0,12 = 0,91 m/sec. 


Taking into account the condition (Vn) < Vi and taking into 
account the friction forces, to determine n using formula 
(18), we will first calculate: 


Q = 2(Virp - Vi, stp) = 2-0, 19 = 0, 38 m/sec; 

R=2 Virp - 3 Vi,5rp = 2:0, 19 = 0, 38 m/sec; 

2Vo + Vs - Q= 2, 32 + 0,50 - 0, 38 = 2, 44 m/sec; 

22 Vo + V+-Q= 2,442 = 5, 8536 m2/sec?; 

4V« = 4-0, 50 = 2, 00 m/sec; 

8V:(Vo - R) = 4, 00(1, 16 - 0, 38) = 4-0, 78 = 3, 12 m?2/sec?; 


Jv, +V. - OY —8V.(V, — R) = 2.733 =1,65 m/sec 





Based on the above calculated calculations, we determine the 
n value: 


_2, 44+1,65 4,09 
2,00 2,00 





= 2,04, n = 2, 04; 


ty =n-ty = 2,04-2,50 = 5,10 sec; from experience ty = 5,10 sec; 
V2,04 = - 0,56 m/sec; 

AH2 = 93, 98(0, 39 + 0, 56) = 89, 28 m; 

from experience AH? = 86,0 M. 

2. According to V experience, we have data: 


Vo = 1,50 m/sec; Ho = 10,0 M; Hrp = 82,5 M; Heaxmax. = 8, 0 M; 
a=545 m/sec; ty = 4, 80 sec; d= 50 mm. 


Next, we calculate: 


V« = 0, 32 m/sec; Dae = (),66 m/sec; = 55,56 sec; 
2d g 


a 

Vip = 0, 66:0,322 = 0, 07 m/sec; V* = Vs - V'rp = 0, 25 m/sec; 
Vi = Vo - V+ = 1,50 - 0,32 = 1,18 m/sec; 
Virp = 0, 66-1,18? = 0, 78 m/sec; 


Vis = Vi - V+ - Virp = 1,18 - 0,32 - 0,78 = 0,08 m/sec; 


Vi, 5rp = 0, 66-0,082 = 0, 00 
V2 = Vi,5 - V+ - Vi, stp = 0, 08 - 0, 32 - 0, 00 = - 0, 24 m/sec; 


Vorp = 0, 66-0,242 = 0, 04 m/sec; 
V2,5 = V2 - V+ + Varp = - 0,247 - 0,32 + 0,04 = 0,52 m/sec; 


V2, stp = 0,66-0,52? = 0,18 m/sec; 
V3 = V2,5- V+ + Vasrp = - 0,52 - 0,32 + 0,18 = - 0,66 m/sec. 


Taking into account }} Vi = 0, we can write: 
Q=2(0, 78 + 0, 00 - 0, 04 - 0, 18) = 2:0, 56 = 1, 12 m/sec; 


R = 2Virp + 3V1, stp - 4Varp - 5V2, stp = 1, 56 - 0, 16 - 0, 90 = 0, 
50 m/sec; 


Vo -R=1,50 - 0,50 = 1, 00 m/sec; 
2Vo + V#-Q=3,00+0,32- 1,12 =2,20 m/sec; 2,202 = 4, 84; 


4V« = 4-0, 32 = 1, 28 m/sec; 8V+(Vo - R) = 2, 56-1, 00 = 2, 56 
M?/sec?; 





J(2v, +V. - Q)° -8V.(V, — R) = 2.28 = 1,51 m/sec 
"*R256L5N) 3,71 
1,28 1,28 





=2,90 n=2,90; 


ty = nt» = 2,90-4,80 = 13,92 sec; no onpity ty = 11,10 sec; 


V2,90 = V25 - 0,40/0,50(V3 - V25) = - 0,52 - 0,11 = - 0,63 
m/sec; 


AH2= 55, 56-0, 88 = 48, 9 m; no onbity AH2 = 49, 0 M. 


4. CONCLUSIONS 

Taking into account the above, in conclusion, we can draw 

conclusions. 

1. The dependence is obtained for calculating the 
maximum pressure during a hydraulic shock in the first 
period of reduced pressure in the event of a rupture of 
the continuity of the flow, taking into account the 
pressure loss due to friction in a long pressure pipeline. 

2. The ratio for determining the number of phases during 
the impact process is obtained. 

3. The obtained ratios are very important for the analysis 
of technical parameters of water hammer. 
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